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We consider the model involving the oscillation of the effective gravitational constant 
that has been put forward in an attempt to reconcile the observed periodicity in the 
galaxy number distribution with the standard cosmological models. This model involves 

> 

00, 

. a highly nonlinear dynamics which we analyze numerically. We carry out a detailed study 

O 

of the bound that nucleosynthesis imposes on this model. The analysis shows that for 

On 

cr 



any assumed value for Q (the total energy density) one can fix the value of f^ar (the 
baryonic energy density) in such a way as to accommodate the observational constraints 
coming from the 4 He primordial abundance. In particular, if we impose the inflationary 
value = 1 the resulting baryonic energy density turns out to be f2b ar ~ 0.021. This 
result lies in the very narrow range 0.016 < f2b ar — 0.026 allowed by the observed values 
of the primordial abundances of the other light elements. The remaining fraction of Q 
corresponds to dark matter represented by a scalar field. 
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I. INTRODUCTION 



The recent observations, in deep pencil beam surveys ][]], showing that the galaxy 
number distribution exhibits a remarkable periodicity of \28h~ 1 Mpc comes as a shocking 
development, since, if taken at face value it would imply that we live in the middle of a 
pattern consisting of concentric two-spheres that mark the maxima of the galaxy number 
density. This, of course, would be catastrophic for our cosmological conceptions. While 
it is true that such periodicity has been observed only in the few directions that have 
been explored so far, it would be a striking coincidence if it turns out that it is absent 
in other directions and we just happen to have chosen to explore the only directions 
in which that phenomenon occurs. Therefore it seems reasonable to assume that the 
periodicity is also present in the deep pencil beam surveys in other directions, thus taking 
us to the concentric spheres scenario. The seriousness of the situation is such that this 
type of scenario has indeed been put forward in a model where the formation of these 
concentric shells is a result of a "spontaneous breakdown of the cosmological principle" via 
a mechanism that results in the appearance of patches filled with the concentric spheres 
pattern, with these patches covering the Universe 0. Again it would still be difficult to 
explain how we did come out living in the center of such a patch (more precisely inside 
the innermost sphere of one such patch). 

The only known way out of this type of scenario is to assume that there is only an 
apparent spatial periodicity that is the result of a true temporal periodicity which shows 
up in our observations of distant points in the Universe and that is mistakenly interpreted 
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as a spatial periodicity ||. The models that have been put forward in order to achieve 
this temporal periodicity involve the oscillation of an effective coupling constant due to 
a contribution to it coming from the spectation value of some scalar field, that actually 
oscillates coherently in cosmic time in the bottom of its effective potential. 

The specific models that have been proposed involve the oscillation of the effective 
electric charge, electron mass, galaxy luminosity or gravitational constant j|,[4||||6|] . From 
these the first two have been shown to conflict with bounds arising from the test of the 
equivalence principle J7|. As for the third scenario, it would seem to involve a large num- 
ber of hypotheses since the galactic luminosity is fixed by the number and type of stars 
present in the galaxy and their respective luminosities, and the latter are themselves 
functions of the standard physics coupling constants that control nuclear reaction rates 
(on the variation of which there are severe limits), and of the transport mechanisms: 
convection, radiation, etc., which are also determined by the standard physics coupling 
constants. Thus it seems that the only way to produce such a model requires the intro- 
duction of "exotic" particles (like axions or massive neutrinos ||) that would act as a new 
transport mechanism and besides that, the hypothesis of a second mechanism that would 
produce the oscillation, an assumption which presumably involves the coupling of these 
hypothetical particles to the hypothetical cosmological scalar field. 

In light of the complicated nature of the alternative scenario, it seems worthwhile to 
carry out a careful analysis of the viability of the oscillating gravitational constant model, 
despite the difficulties that seem to appear when confronting the predictions of the model 
with other experimental data. We will address these difficulties below. 
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The oscillating G model is based on a cosmological, massive scalar free field, that is 
nonminimally coupled to curvature and whose oscillations in cosmic time result in the 
oscillation of the effective gravitational constant. The model is governed by a system 
of nonlinear, ordinary, differential equations that we shall integrate numerically without 
analytical simplifications. In previous related works P,0fj], only approximate solutions 
have been explored as in the work of Morikawa 0] . In this work, the dynamical equation 
for the scalar field is linearized into a Bessel equation. From this, the Hubble parameter 
for a flat matterless Universe is evaluated. The approximations result in a precision of 
2% in the Hamiltonian constraint. While such procedure is an extremely simple way to 
study the implementation of the model it is however not totally satisfying. In addition 
to the analytical simplifications Morikawa made, we emphasize that his "matterless" (no 
baryonic or radiation energy density) Universe assumption is not useful for the study of 
the constraints we discuss below, notably the one imposed by nucleosynthesis. As we shall 
analyze in Sec. IV, it is actually the presence of matter that makes it possible to satisfy 
such a constraint. Furthermore, it is obvious that Morikawa's approach is no longer valid 
in early times where the linear approximations break down. Hence, this method cannot 
be used to describe the behavior of cosmological variables during the nonoscillating phase. 

The data arising from the Viking radar echo experiments in addition to those related 
with the limits on the Brans-Dicke parameter || impose certain bounds on the value of 
G/G. There is an apparent conflict between these bounds and the value they should have 
for reproducing the galaxy-amplitude counting. However, this problem can be overcome 
by imposing the unnatural extra requirement that our Galaxy is at a "fortunate phase" 
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H where the scalar field is swinging very close to zero at our particular place and time 
in the Universe. There are additional bounds on the oscillating G model arising from Big 
Bang nucleosynthesis, but we will argue that these cannot be adequately addressed in 
the fashion described by Accetta et al. || . This constraint is obtained by considering the 
Helium abundance limits which depend on the neutron-proton ratio at the temperature 
of freeze-out. This temperature is itself obtained from the condition that the Hubble 
parameter equals the nucleon to proton weak reaction rate. 

Crittenden and Steinhardt B, based on a previous analysis by Accetta and collabora- 
tors that resulted in the bound AG/G < 0.4 ||, argue that the nucleosynthesis constraint 
is so stringent that it practically rules out the oscillating G model unless we assume a 
"fine-tuning" within the oscillations of the scalar field. 

The problem with employing the analysis of Accetta et al. [f| directly to this model 
is that they study the change ocurring in the Hubble parameter as a function of the 
temperature when one changes the value of G, but fails to take into account the fact that 
the model implies an equation for the Hubble parameter that contains terms other than 
the simple ones used by the authors. These terms are associated with the contributions 
from the scalar field to the effective energy density. 

In this paper we carry out an analysis of the nucleosynthesis bound taking into ac- 
count these extra terms and demonstrate that indeed such a "fine-tuning" is needed and 
moreover that it is possible. However, we will argue that this "fine-tuning" is not of the 
kind that should result in the dismissal of the model, but rather a natural adjustement of 
the initial conditions that will lead to the observational data extracted from our Universe 
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today. In other words, scientific models that require a very precise choice of the numerical 
value of the initial conditions in order to reproduce a given qualitative behavior of the 
observational data, are models that would be consider unnatural and the choice of the 
specific initial data is justifiably described as "fine-tuning." However, models that require 
a very precise choice of the numerical value of the initial conditions in order to reproduce 
a specific numerical observational data cannot be considered as unnatural, especially if 
for every conceivable value of the observational data (at least in some range) there is a 
corresponding value of the initial data. In this type of models, the particular "preferred" 
value at the initial data is just the result of a 1 to 1 correspondence between initial con- 
ditions and final outcome. We will argue that the present model could be of the latter 
type. 

The organization of the paper is as follows: In Sec. II the model is described and the 
basic equations of evolution for the fields and matter are derived, in Sec. Ill the numerical 
implementation is discussed together with the error estimation analysis, in Sec. IV the 
results of the numerical integrations are analyzed and finally in Sec. V we give a brief 
discussion of the main features exhibited by the model, their physical significance and the 
overall viability of the model. 

II. FORMULATION OF THE MODEL 

We will consider a model in which the effective gravitational constant becomes depen- 
dent on cosmic time due to a contribution to it coming from the spectation value of a 
scalar field. This can be achieved by considering a scalar field nonminimally coupled to 
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gravity. One of the simplest models of this kind is obtained by taking a Lagrangian as 



follows 



C 



i(v0) 2 + y(0) 



mat 



Here Go is the Newton's gravitational constant; £ stands for the nonminimally coupling 
constant; and V(4>) is a scalar potential to be specified later (see Sec. III). In this model 
we are also including an schematic matter Lagrangian £ mat . Equation ([]]) shows that the 
introduction of the coupling term is equivalent to considering an effective gravitational 



constant which explicitly depends on the scalar field: 



G, 



cff 



1 + 16ttGo^ 2 



(2) 



The gravitational field equations following from the Lagrangian ([I]) can be written as 



where 
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(3) 



mat ) i 



Tj£ = G cS (4£Tf + T s f + TZ 



;(V0) 2 + V(0) 



(4) 
(5) 
(6) 



The energy-momentum tensor of "matter" T^t wm De composed of a combination of two 
noninteracting perfect fluids, one corresponding to pure baryonic matter (z = 1) and the 
other one representing a pure radiation field (i — 2): 



J mat J bar ' f / / 

i=l,2 



[(pi + ejUW+pig 



(7) 



which posseses the symmetries of the spacetime. The scalar field will also be assumed to 
posses these symmetries. 

Finally, the equation of motion for the scalar field becomes 



□0 + 2£<pR 



dV{(p) 



We will focus on the Friedmann- Robert son- Walker (FRW) spacetimes which describe 
isotropic and homogeneous cosmological models 



ds 2 = -dt 2 + a 2 (t) 



dr 2 



1 — kr 2 



+ r 2 d9 2 + r 2 sin 2 



(9) 



where a(t) is the scale factor and k — 1, 0, — 1. 

Our purpose is to study the behavior of the solutions of the gravitational, matter, 
and scalar field equations for the FRW line element (|9|). Since these equations are highly 
nonlinear, it is a difficult task to find analytic solutions; therefore, we will approach the 
problem via a numerical analysis. 

One of the equations we find is the Hamiltonian constraint 



d 2 k 8 _, 



(10) 



The dynamical equation for the single gravitational degree of freedom is 



a a 2 k ^ ( 1 . 

- + 2- + 2- = 47rG Q (E--S) , 
a cr cr \ 6 
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where 



E 



G 



■ a 



e+-0 2 + V(0)-12e#- 
G L 2 a 



(12) 



is the total effective energy density and 
3G e ff 



S 



p + U 2 - V(cP) + 4£ U 2 - #- - 0D( 
2 V a 



(13) 



These source terms contain contributions from the three parts of the total energy- 
momentum tensor given in eqs.(^)-(0). 

Finally, the equation for the scalar field @ can be written explicitly as 

+ 30- + = 16nG ^(E - S) , (14) 

where we have replaced the scalar curvature in terms of the energy-momentum tensor 
quantities E and 5". 

It is worth mentioning that the contribution of the scalar field to the expression T^. v = 
[see Eq. (^)] vanishes identically and thus the energy-momentum of the ordinary matter 
satisfies T^ t . v = 0. We will moreover assume that the two perfect fluid components 
(baryons and photons) do not interact among themselves, thus each of their corresponding 
energy-momentum tensors is separetely conserved leading to 

e i + 3(e i + Pi )- = . (15) 
a 

Equation (|15"D integrates immediately with respect the scale factor like in the standard 
cosmology case. We find then 

e = e bar + e 7 = ci (^j + c 2 , (16) 

P = Pbar + P 7 = y (^j with a := a(t = t ) . (17) 

Here we have assumed an equation of state p 1 = e 7 /3 for the radiation part, whereas 
Pbar = for the corresponding baryonic component. The first term of (0) represents thus 
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the pure baryon energy density, while the second one the radiation contribution alone. 
The constants c\ and C2 are fixed by the "initial" conditions (i.e., the conditions today). 
In particular (at t = to) we will be assuming a baryonic energy density of one tenth of 
the critical value (at least for the first numerical experiment, but latter we will use C\ 
as an adjustment parameter) and the radiation energy density corresponding to the 2.73 
K cosmic background radiation (CBR) (see also Sec. III). Hereafter we shall refer as to 
matter the combination of both fluids. 

With the aim of reducing the field equations into an initial value problem consisting 
of a system of first-order differential equations, we shall rearrange these conveniently. 
Here we present the final form of equations with source terms containing no second-order 
derivatives and introduce better-suited variables: 

P Q = -a(t) = ~, P<P = <P, (18) 



where 

~a(t) 



a(t) := In 

The dynamic equation ([□]) then takes the form 



(19) 



P a -P a 2 = ^nG (E + S) , (20) 



where we have used the Hamiltonian constraint (10) in order to eliminate from eq.flTTD 
the term proportional to k. Notice that P a = —H(t). 

Introducing eqs.(^) into the scalar field equation we obtain 

P - 3P a P^ + = 167T G ^(E - S) . (21) 
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The source terms then take the following form: 



E 



G 



cff 



Go 



e + -P/ + V{4>) + 12^P a 



(22) 



and 



S 



G 



cff 



1 + 1927rG' eff ^ 2 2 G c 
+ 4£ (<pP^P a + P<f> 2 -< 



V+ -P^-V{4>) 

d ^jf) + ^G e^E 



(23) 



In obtaining the source fl23|) from (|T3|) we have used the scalar field equation (pl|) in order 
to eliminate the term with □</>. 

Further analysis of this model requires the numerical integration of the field equations 
under appropriate initial conditions. This will be performed in the following sections. In 



the Appendix we provide the dimensionless form of the above equations. 



III. INITIAL CONDITIONS AND FIXING OF PARAMETERS 



The choice of our variables fixes in advance the values of a and P a at present time. 
Then, the initial conditions are 

a\t=t = , (24) 

P a \t=t = -l ■ (25) 

It is worth emphasizing that by initial conditions we will mean throughout the paper the 
value of the field variables at present time and not their corresponding value near the Big 
Bang. 

11 



As stated before, the initial condition for e corresponds to the value of the (baryon 
plus radiation) energy densities today and then: 

e := e\ t=t0 = e£ ar + e J = a + c 2 . (26) 

We fix C2 by the value T = 2.73 K from the CBR and choose for the moment a value of 
ci corresponding to f2bar = 0.1. Thus, the values of the constants appearing in Eqs. ([TB]) 
and ( |T7| ) are 

ci = 0.1e c , (27) 

3c 2 if 2 

c 2 ~ 4.2 x l(T 5 e c with e c := — . (28) 

87rGo 

Here e c stands for the critial energy density in terms of the current value of the expansion 
rate H . 

This choice will leave Q (see below) as a free parameter which is adjusted in order to 
determine the initial conditions on which result in a model respecting the observed 
4 He abundances. As we mentioned before, we will be able to follow (Sec. IV) a different 
strategy that consists in imposing ft — 1 and using c\ (i.e., fi bar ) as an adjustement 
parameter. 

The Viking data experiments constraint Gj (G-ff)|toda y to be less than 0.3/i _1 0. This 
imposes a bound on <p. The most conservative and the one we choose is 



|t=*o 



. (29) 



However is straightforward to explore a less conservative initial condition. 

The values of the initial scalar-field's amplitude and coupling constant are obtained 
from the observational data. The observed redshift-galaxy-count amplitude (see HH) and 
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the Hamiltonian constraint (|TI]) at t provide two algebraic equations which determine 
the values of 0o and £ once we choose a value for fl and f2bar- When assuming a harmonic 
scalar potential V(<f)) = m 2 (f) 2 , these equations are 

A = -gf^ , (30) 

where eo is the initial matter energy density in units of e c ; Q := Eo/e c and u; 2 oc m 2 is 
the oscillation frequency in units of Hq 1 . 
Solving for O an d £ we obtain: 

= Am 2 f 32) 

s 167r[^ e +a;(0-eo)] ' v ' 

^ = 16^ [u- M ■ (33) 
The model will explain the observed galaxy distribution if ^lo ^ 0.5 0. The remaining 
parameter to be adjusted is u which is obtained from the observed galaxy periodicity of 
128 Mpc hr 1 . 

In order to test our numerical code we have first restricted ourselves to the standard 
cosmology where known analytical solutions exist. The field equations have been solved 
by means of a fourth-order scheme with an adaptive stepsize control. We have performed 
the integration of the equations with respect to time and also with respect to the variable 
a. The time integration produces relative errors on the dynamical variables, like the scale 
factor, energy density, and Hubble parameter which are of the order of 10~ 8 . 

The choice of the variable a as integration parameter allows us to explore the evolution 
at very early stages (a < region) while keeping the relative errors small. Indeed we stop 
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the integration at a ~ —25, a bit beyond the value at which the nucleosynthesis takes 
place, the calculations having started at a = 0. 

While the integration of equations for flat (k = 0) and hyperbolic (k = —1) Universes 
can be performed for an arbitrarly large value of a, for a closed Universe (k = 1) the inte- 
gration makes only sense for a < ct max . The limiting value corresponds to the maximum 
size reached by the Universe and beyond which it starts recollapsing. We mention that 
the regions a > a max are indeed not very interesting, first because the physics associated 
with them can be inferred from the a < a max branch and then because no observational 
bounds arise from that region. 

We use the Hamiltonian constraint as a test on the accuracy of the procedure. It is 
applied at every integration step and implemented by defining the deviation parameter 
(see the Appendix) as, 

_ (H/H f - (g - 1) e- 2 " 

For an infinite accurate integration this parameter would equal one and the deviation 
from this value indicates the degree to which Eq. fllTf ) fails to be satisfied. 

Figure p] shows the correlation between the relative errors found for the cosmological 
variables and that for the Hamiltonian constraint (i.e., in A) as a function of time for the 
Einstein-de Sitter Universe. The deepest peaks (infinite precision) indicate the location 
of initial data. Only for convenience, in such points the precision has been arbitrarily set 
to be ~ 10 16 . 

We have also verified this type of correlation in closed and hyperbolic Universes. 
Because, no analytical solutions are known for the oscillating models, we use Eq. (|34|) 
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systematically (the "internal test" ) in order to verify the accuracy of the results. 



IV. NUMERICAL RESULTS 



As we will see below, for any given value of f2 it is possible to fix f2bar such that the 
correct value of the He abundance is recovered. Figure shows the behavior of in some 
range of past and future epochs with Q = 1 and the corresponding ~ 0.021. As the 
scale factor increases, the amplitude is damped due to redshifts while the frequency of 
oscillation (with respect to a), which is inversely proportional to H 2 ||, grows (see Fig. 
|3] for the behavior of the Hubble parameter). The amplitude and oscillation frequency 
will reach maximum and minimum values respectively as the scale factor decreases before 
entering into a stage (discussed below) in which the scalar field almost vanishes. This 
includes a late era in which matter is dominant over the scalar field energy density and 
an earlier one in which this is the opposite (see Fig. |j). 

We introduce effective energy densities by 

E ( := U^P^ , (35) 

(jr 







\P, 2 + mV 



(36) 



^mat := ^r-e . (37) 

Figure f| shows how these fractions of the total effective energy density (depicted in Fig. 
H) vary with a. 

We emphasize that E^ (a "coupling energy density" ) is not positive definite and since 



it contributes to the total effective energy density E [see Eq. ( plq) ], there are some regions 
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where the fraction E^/E is negative (dash-dotted line) while E^/E and E mat /E exceed 
one (dashed and solid lines, respectively). In particular, when comparing Figs. |2| and |] 
we note that when the amplitude \4>\ is large (in the region a G [—4, —3]), the net scalar- 
field's contribution to E represented by + is small because the negative "energy 
density" E$ compensates the contribution from E^. Indeed, it is the matter contribution 
E msA , which becomes dominant in this era. 

When the scale factor is still small (a G [—8,-5]), \<f>\ falls down dramatically by 
entering what we called the "fine-tuning era" (see Sec. IV A) where E mat dominates 
completely. The effective gravitational constant which was reduced in ~ 38% of its current 
value Go during the maxima of \<p\, recovers its normal value again (Fig||). Finally, for still 
smaller a, the scalar field becomes dominant again, resulting in G c k — > as we approach 
the Big Bang singularity. The Hubble parameter decreases monotonically in the early era 
(a G (-oo,3]) (Fig. §. 

A. The nucleosynthesis bound and "fine-tuning" 

In previous investigations, Crittenden and Steinhardt || have brought attention to a 
couple of much more severe constraints on the oscillating models that those considered 
before by Morikawa j|] and Hill et al. [[| . The first of these constraints regarded the bound 
imposed by the Brans-Dicke parameter tests which are even more stringent that the one 
coming from the Viking radar echo experiments. However, the authors argue that it can 
be eluded by assuming a "fortunate phase" for the oscillation of the scalar field when the 
nonminimally coupling function varies quadratically with the changes of 0. The second 
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new bound they mentioned arises from the fact that prior to the onset of its oscillatory 
behavior, the amplitude of had to be small enough in order to prevent \AG\/G > 0.4 
at nucleosynthesis as described by Accetta et al. ||, but at the same time, the evolution 
of cf) has to be such that its present amplitude is large enough to accomodate the value 
of Aq needed to generate the observed peaks in the deep pencil survey Q. Crittenden 
and Steinhardt suggested that both conditions might be possible by "fine-tuning" the 
amplitude of <fi at nucleosynthesis. This amplitude might then be amplified by the effect 
of the curvature during the matter-dominated era and afterwards damped by the redshift's 
effects. 

Although the nucleosynthesis bound can be thought as imposing a limit on the vari- 
ation of G, this is not enterely precise. The true bound that comes from nucleosynthesis 
is a limit on the deviation of the expansion rate of the Universe from its value given by 
the standard cosmological model. There exists a very narrow region for which the expan- 
sion rate of the Universe and the transition rate for the weak interaction, which convert 
neutrons to protons, traduce themselves into a freezeout temperature that reproduces the 
observed 4 He abundance (see [TO] for a review). In the standard cosmology framework, 



for a radiation dominated era, the ratio between those rates, is 

kT 



(3* 



H V0.7 MeV, 

Thus the temperature at which those weak interactions freeze-out the n/p ratio is nTp ~ 
0.7 MeV. At this temperature, the neutron-proton-ratio is given by its equilibrium value 



n 



exp {-Q/T F ) « 1/6 , (39) 
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where Q is the mass difference of neutrons and protons. This ratio can decrease to 1/7 
if we take into account the "natural" neutron decay due to weak interactions. This value 
predicts approximately one 4 He nuclei for each 16H nuclei. Therefore the primordial 
abundance of 4 He-H ratio gives Yp ~ 25%, i.e., very close to the observed value, varying 
in some small amount depending on the number of neutrino species considered. 

Let us emphasize that one of the greatest achievements of the standard cosmology is 
the predicted light element abundance that we observe today. Therefore, the oscillating 
model would be considered as viable if it preserves also these nice features. 

Because of the exponential relation Q55D, a small variation of the freezeout temperature 
will produce a large deviation of the observed neutron-proton ratio. Deviations of the 
freezeout temperature will arise as a result of a small change of the expansion rate. The 
greater the expansion rate, the greater the freezeout temperature and so larger the 4 He 
abundance. This is the main reason why alternative cosmological models can fail. 

We showed in the previous section that the initial condition O an d the value of £ 
depends upon Ao , 00, Q, and f2bar- The parameters to and Ao are fixed in order to satisfy 
the observed galaxy periodicity and the galaxy peak amplitude whereas <fi = prevents 
conflicts with the Viking experiments constraint. So we can try to adjust Q or Q^ aT in 
order to ensure the satisfaction of the nucleosynthesis bound. Initially we considered 
fixing fibar — 0.1 and then adjusting Q. Another perhaps more natural possibility (see 
discussion) is to take a value of Q (in particular Q = 1 as suggested by inflation) and then 
adjusting fibar- As we will see below, we have done this obtaining f2 b ar ~ 0.021. 

In the process of integrating towards small values of a we find initially the following 
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behavior: for some choices of f2 bar , the value of goes monotonically to +00, or to —00 
for some others. This suggested to us that there exists some value of f2bar at which the 
transition from one behavior to the other takes place. We found that for Q = 1 this 
"fine-tuned" value turns out to be £\ar ~ 0.021. Figure |7] shows the behavior of <f> for 
three values of fl^ aT about this point. 

This transition point, traduces in a very special "initial" conditions for which the 
amplitude is "squeezed" to zero, before the onset of the oscillatory behavior. The 
search of that transition value is what we call "fine-tuning. " Our numerical experience 
indicates that it is possible to generate a finite region during which <fi ~ before growing 
or decreasing monotonically (in the direction of decreasing a) and that the extent of that 
region depends on the improvement of the "fine-tuning" . 

Different "experiments" in "fine-tunings" are depicted in Fig. [g. The larger the 
plateau for which G c r — > Go (i.e., for which — > 0) the closer the freezeout temperature 
is to the standard cosmology prediction 0.7 MeV. We have found that fl^ = 0.021 + e 
results in a freezeout temperature which agrees with 0.7 MeV and therefore giving a 4 He 
abundance that approximates best the observed value (see Fig. §). 

In light of the extreme sensitivity of Tp to the value of f2bar we are not able to improve 
the "fine-tuning" due to the limitations in the numerical precision. However, even with a 
noninfinitely precise "fine-tuning" the model is able to recover the 4 He abundance from 
nucleosynthesis. Figure [5] shows the freezeout temperature (~ 0.7 MeV) predicted by the 
best "fine-tuning" we explored. 

The Q = 1 scenario (with £\ar ~ 0.021) corresponds to an age of the Universe of 



19 



~ OSHq 1 (see Fig.[10|) . Figure [□] shows a typical curve of the Hubble parameter as a 
function of redshifts. The transition between the matter and radiation epochs is clearly 
appreciated in Fig. |T2] [see also Eq. fllfip]. 

Only for completness we mention that, for example, if we choose £\ar = 0.1 and 



fine-tune this results in Q ~ 1.67 which corresponds to a closed Universe. Figure [13 



shows that the age of the Universe based in this scenario is about 0.7 Hq 1 . Somewhat 
younger than the age predicted by a Universe that includes only baryonic matter and 
radiation energy densities (a hyperbolic-standard-cosmology model) which corresponds 
to ~ 0.9HQ 1 . Figure [H] shows the freezeout temperature (~ 0.8 MeV) predicted by the 
best "fine-tuning" we explored in the fibar = 0.1, Q ~ 1.67 scenario. This temperature is 
very close to the one predicted by the standard cosmology (~ 0.7 MeV). This examplifies 
the fact that for any choice of Q (or fi bar ) the freezeout temperature of standard cosmology 
can be recovered by adjusting the other parameter f2bar (or Q). 



V. CONCLUSION AND DISCUSSION 

At first sight the model could not be less appealing: a "fortunate phase" and an 
incredibly precise "fine-tuning." However, we must judge it in the proper context by 
contemplating the alternatives. 

First the observed periodicity might be a statistical fluke, then of course that would 
be the end of the story. However as we argued in the Introduction it would be quite a 
coincidence that we just happened to look at the couple of directions that exhibit such 
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patterns. In any event this is a matter that only further observations will answer. 

Next we have the "spontaneous breaking of the cosmological principle" O, clearly a 
major departure form our cosmological conceptions which besides that, requires something 
even more fortunate than the "fortunate phase" of the oscillating G models: We happened 
to be "fortunate" enough to be born in a galaxy which happens to lay in the middle of a 
concentric collection of shells of maximal galaxy density. 

And finally the galactic luminosity oscillation model, which as we said at the beginning, 
more than a model is a vaguely specified scenario which nevertheless seems to require 
at least two new hypotheses: a new type of star cooling mechanism, and also (as the 
other alternatives do) a driving oscillating cosmological scalar field to turn on and off 
that mechanism periodically in cosmic time. Needless to say that once the scenario is 
implemented with a specific model, unforsen new bounds might also have to be overcome. 

In this light, the oscillating G model does not look as clearly dismissable. Furthermore, 
while it is true that the fortunate phase will have to remain such, for the other problem, 
"the fine-tuning," we will argue below that there are scenarios in which this problem is 
not present. 

To start, we must stress that while the "fine-tuning" is completely unnatural when 
approached, as we have, from the present to the past, when looked from the opposite 
and more natural direction, the situation is quite different. In fact all that seems to be 
required is for some mechanism to drive the scalar field to an extremely low value before 
the era of nucleosynthesis, for then our calculations show the field will remain at that 
value up to and beyond that era so will have G e s ~ Go and then the success of "Big Bang 
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nucleosynthesis" will be recovered naturally. The amplitude \<j>\ will later be amplified by 
the curvature coupling precisely before the onset of the oscillatory behavior. 

Thus it is possible that starting from an arbitrary value of near the Big Bang, a 
mechanism (no different from that required to solve the other problems of the standard 
Big Bang model, including inflation itself) would drive the scalar field to a value near 
zero, at which it will remain until just before H ps m when the amplification and then 
oscillations would occur. This type of scenario might be combined with the inflationary 
prediction Q = 1 and a corresponding "fine-tuning" on fibar- This has been done and 
the "fine-tuning" yielded a value f2bar ~ 0.021, surprisingly in the very narrow range 
0.016 < f2 bar < 0.026 that results in a successful nucleosynthesis of the light elements 
other than 4 He. But the point is that this would really be no "fine-tuning" at all (if 
looked in the right perspective), because all that it will mean is that, given the physical 
constant precise values, the inflationary mechanism will ensure that the energy densities 
of the various components, scalar field and ordinary matter, add up to Q = 1, which will 
correspond then to a Universe at our time with precise values of the densities, expansion 
rate, etc. In particular the precise current value of the scalar field amplitude and phase 
arises from a particular precise value of the parameters at early times, among them the 
baryon content of the Universe. The fact that the corresponding baryonic density today 
turns out to lay in a very narrow range consistent with the light element's nucleosynthesis 
suggests that the model should be taken seriously. In this respect we would like to point 
out that the a priori probability for this happening just by chance is about 1 in 100 (the 
range [0.016, 0.026] which is observationally allowed for f2 ba r represents approximately 1 
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part in 100 inside the numerical range allowed in principle [0, 1]). 

In view of the previous arguments we may conclude by saying that the "fine-tuning" 
should be seen as simply recovering the "initial" conditions corresponding to the observa- 
tional data and that, moreover, this has resulted in the specific prediction f2b ar ~ 0.021. 
Therefore, this is certainly the most atractive of all the models considered in order to ex- 
plain the observed periodicity in the galactic distribution, and should also be considered 
as a missing mass model with the scalar field playing the role of dark matter which is, 
however, indirectly observable in the oscillation of the galactic distribution. 

APPENDIX: DIMENSIONLESS FORM OF FIELD EQUATIONS 



1. Einstein's dynamical equation 



It is easy to check that when restoring factors of c and introducing characteristic 



lengths, Eq. ( |20"D becomes 



P n = Pi 



I \e + s 



(Al) 



with 



t := tH n 



P a :~- 



E :-- 
S :-- 
e r :-- 
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(A2) 
(A3) 
(A4) 
(A5) 
(A6) 
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Here a dot over a tilde means derivation with respect the dimensionless time t. 



2. Scalar field's equation of motion 



Using fl£3 



and introducing 



Pa :-- 



(A7) 



Eq. ( pi]) can be written 



= 3P ot P <i ,-^-V' (</>)- 6^ <f> 



5-P 



(AS 



The scalar field now turns to be a dimensionless quantity and 



(A9) 



Thus, in the case of a scalar potential V(<f>) = A0 n , A has units of energy density. In 
particular if A := m? and n = 2, the dimensionless harmonic frequency of oscillation is 
given by 

0,:=- = ^. (AID) 
Moreover, with such a choice of units, the coupling constant £ is also dimensionless. 



Finally, the dimensionles source terms read explicitly 

3G 



S 



P+ T P0 -W) 



1 + 192ttG£ 2 <p 2 



+ 1 + 1927rG£ 2 </> 2 
P = G 



47T ~ 2 ~ ~ 
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1 + 167r{. 



(All) 
(A12) 
(A13) 
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Equations (TATp , flAg), (|A7p, ( gg ) with source terms QCT])- (pnB|) and the equation 
of conservation of energy ( |A20| ) are then the complete set of equations best suited to 
be solved numerically. With a trivial manipulation we write these equations in terms of 
derivatives with respect to a instead of t. 

3. The Hamiltonian constraint and the deviation parameter 

With these notations, the dimensionless form of the Hamiltonian constraint ( JlQ| ) reads 



p2 + 



a a 2 m 



E 



(A14) 



Now, at t = to this becomes 



1 + 



n 2 H 2 



En . 



(A15) 



We can now replace k in Eq. ( |A14j ) in terms of hereabove initial conditions to get 



P*+(E -l)e- 2a -E = , 



(A16) 



where we have employed our variable a instead of a. 

We can introduce also the dimensionless deceleration parameter q(t) in terms of sources 



m 



i 



in 



E + S 



(A17) 



At t = in, we can also rewrite (|A17|) in the form 



Eft = 2q — Sq 



(A18) 



Moreover, the deviation parameter introduced in (|34| ) take the form 
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A :-- 



E 



-2a 



(A19) 



where we stress P a = —H(t)/H and fl = E . 



4. Conservation equations of matter and radiating fields 



Equation (|T^) becomes 



li = 3(ei+pi)P Q 



(A20) 
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FIGURES 

FIG. 1. Relative errors as a function of dimensionless time: solid 
lines(superposed): log 10 | (a(t) an - a(i) num ) /a(*) an |, log 10 | (d(t) an - d(t) num ) /d(t) an |; dashed 
line: log 10 | (a(t) an - a(t) num ) /a(t) an |; dash-dotted line: log 10 | (P Q (t) an - P a (t) num ) /P™\; dot- 
ted line: log 10 | ((ea 3 ) an — (ea 3 ) num ) /(ea 3 ) an |; dot-dashed line (bottom): log 10 |l — A|. The initial 
conditions located at to := O-^Hq 1 ■ 

FIG. 2. Fine-tuned scalar field amplitude as a function of ln[a/ao] for a flat Universe (f2 = 1) 
with r2b ar = 0.021 012 641 182 345 and Ao = 0.5 at the onset of oscillations. At present time 
(q = 0) the initial amplitude is 0o ~ 3.288 x 10~ 3 and £ ~ 6.267. Computations were stopped 
at a ~ 1.5 . 

FIG. 3. The Hubble parameter obtained by using the same initial values as in Fig. | and 
H(a = 0) = H . 

FIG. 4. Starting from the same initial values as in Fig. §, this fi gure depicts the fractions 
of effective energy densities: E ma _ t (solid line), (dashed line) and, E^ (dash-dotted line) (see 
text for definitions). 

FIG. 5. The evolution of the total effective energy density resulting from the numerical 
integration of the field equations of Sec. II (see also the Appendix), with initial values as in Fig. 

1- 

FIG. 6. Behavior of the effective gravitational "constant" in units of Newton's constant Gq. 
The initial values have been chosen as in Fig.Q. 
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FIG. 7. Behavior of the scalar field for three different values of f^bar within the = 1 scenario. 
The solid line represents the best "fine-tuning" obtained with Qbar = 0.021 012 641 182 345 (see 
Figj2|); the dashed line for Qbar = 0.022, and the dash-dotted line with Qbar = 0.020. 

FIG. 8. Effective gravitational "constant" in units of Newton's constant Go for differ- 
ent fine-tunings: solid line obtained with f^bar = 0.021 012 641 182 345, dashed line for 
ftbar = 0.021 012 641 18 dash-dotted line with ft bar = 0.021 012 6. 

FIG. 9. Expansion (solid line) and transition (dashed line) rates in terms of the blackbody 
temperature. The asterisk depicts the freezeout temperature ~ 0.7 MeV at which nucleosyn- 
thesis takes place as predicted by the standard cosmological models. The cross-point of both 
curves indicates the corresponding freezeout temperature ~ 0.7 MeV for the oscillating model 
of previous figures. 

FIG. 10. Scale factor of the oscillating model of previous figures in units of its value today as 
a function of cosmic time (in units of H^ 1 ). The present time to has been placed at t = IHq 1 . 

FIG. 11. The Hubble parameter as a function of redshifts. The initial conditions of Figs. || 



and 1C determine the initial value H/Hq = 1 at z(to) = 0. 



FIG. 12. The combined matter-radiation (solid line) and pure radiation (dashed line) energy 
densities with initial values as in Fig. ^. 

FIG. 13. Scale factor for the closed oscillating Universe with 0, ~ 1.67 and f^bar = 0.1 (solid 
line). The dashed line corresponds to the case with no scalar field (standard cosmology). As in 
Fig.lC, the present time to has been placed at t = IHq . 
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FIG. 14. Similar to Fig.|| for the closed Universe with Obar = 0.1 and Q ~ 1.67. Here the 
predicted freezeout temperature closely agrees with ~ 0.7 MeV. 
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